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ABSTRACT 


This study presents the results of the aerodynamic 
design of a transonic axial—flow fan for an aircraft air-— 
conditioning unit. The design is based on fan design=—point 
specifications and the experimental test results for a 
model of a proposed fan-inlet. 

The fan-inlet model was tested under flow conditions 
that gave Reynolds number similarity to the design inlet 
conditions. Three—dimensional pressure probes were used to 
measure the model discharge flow conditions at eight sta- 
tions. The flow was found to lack axial symmetry. There 
fore, several internal modifications were made. Complete 
axial symmetry was not achieved for a reasonable loss in 
total pressure. It was necessary to circumferentially aver= 
age the velocity ates wien and mass—average the total~ 
pressure loss in order to provide a basis for the fan 
design. 

Concurrently another design based on uniform flow and 


model losses was made. The results were nearly identical. 





Thesis by: Ronald William Pyle and Milton Ward York entitled 
limeansonic Fan For Aircraft Air-Conditioning Unit. 
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Ll. Introduction 

An axial-flow fan may be distinguished from an axial- 
flow compressor in that a fan is generally classified as a 
high capacity, low head, single-stage unit which operates 
at a high, constant speed. Although used for different pur— 
poses than a compressor, an axial-flow fan may be designed 
or analyzed with the available theoretical and experimental 
methods for axial-flow compressors. 

One of the principal assumptions made in the aerodynamic 
design of axial-flow compressors is that of axially symmet- 
ric flow. This condition is primarily assumed to reduce the 
general flow equations from a complex three=—dimensional to a 
simplified two-dimensional system. At a given radial loca- 
tion the physical requirement of this assumption is for 
circumferentially constant flow conditions ahead of any 
blade row. 

The single-stage transonic fan, with its relatively 
higher flow capacity and pressure ratio over conventional 
subsonic fans, has been successfully used in modern air- 
craft systems. However, for compact systems the transonic 
turbomachine becomes more demanding of inlet flow uniform 
ity. Thus, proper matching of an inlet=<scroll to a fan 
stage becomes more critical for a given set of dimensional 


restrictions. 


The purpose of this paper was the aerodynamic design 
Of a single-stage transonic fan for an aircraft air- 
conditioning unit. The design was to be based on a list of 
design=point specifications and also on a proposed fan in-= 
let that complied with the geometrical design envelope. 

A double=-scale model of the proposed fan inlet was tes-= 
ted. Internal modifications were made in an attempt to 
attain uniform, axially symmetric flow at a location just 
ahead of the required rotor position. Final model test 
results were used to predict the losses through the fan in- 
let and the flow conditions ahead of the fan rotor. The 
aerodynamic fan design was based on these conditions. 

A similar fan design was made-which was based on assumed 
uniform, ascii? symmetric flow ahead of the rotor. Thus, 

a comparison could be made between the two designs, and 
also with the design specifications. S 

The authors wish to express their appreciation for the 

guidance and assistance given them by Professor M. H. Vavra 


of the Department of Aeronautics of the Naval Postgraduate 


School. 
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List of Symbols 


A 


= 


s 


a 


e 


> & 
wo 


area, fers 

number of blades 

blade height, in 

specific heat, BTU/1b,, wei 
blade chord, in 

diffusion factor 

driving power 

Reynolds number factor 
area multiplier 
gravitational constant, 32.2 £t/sec* 
total enthalpy, BTU/I1b, 
incidence angle, degrees 


incidence angle of uncambered blade section, 
degrees 


conversion factor, 778.2 ft lb/BTU 

weight flow blockage factor 

correction factor in incidence-angle relation 
correction factor in deviation=angle relation 
characteristic length, in 

Mach number 

factor in deviation—angle relation at OT = 1 
revolutions per minute, RPM 

Reynolds number 


slope factor in incidence—angle relation 


(x DW RS 


am 


polytropic exponent 

total pressure 

static pressure 

dynamic pressure 

theoretical average dynamic pressure 
gas constant, 53.34 £t 1lb/1b,, Rn 
radius, in 
radius, inner and outer, in 
entropy, BTU/1b,, R 

blade spacing, in 

total temperature, °R 

static temperature, oe 

blade maximum—thickness—to—chord ratio 
peripheral speed, ft/sec 

velocity, £t/sec 

flow rate, lby/sec 

equivalent flow rate, isi? 

pressure ratio parameter 

distance from annulus outer wall, in 
coordinate along axis, in 


discharge coefficient 


air angle, angle between air velocity and axial 
direction, degrees 


ratio of specific heats 
deviation angle, degrees 


density coefficient 
inlet total-pressure-loss coefficient 


LO 


mes Sarirrges? 


polytropic loss coefficient 

adiabatic efficiency 

pitch angle in model discharge annulus, degrees 
viscosity, lb,/sec ft 

density, slugs/£t> 

solidity, ratio of chord to spacing 
dimensionless flow function 

blade camber angle, degrees 

yaw angle in model discharge annulus, degrees 
angular velocity, radivans/sec 


fan total-pressure-loss coefficient 


Subscripts: 
0 station ahead of inlet 
au station at retor inlet 
2 station at stator inlet 
2=—D low speed two-dimensional cascade 
3 station at stator exit 
10 blade maximum—thickness—to=—chord ratio of 10 per= 
cent 
ad adiabatic 
av average 
ax axial 
c compressor 
design 


model discharge 
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H horizontal 

h hub 

I model entrance 

j summation index 

M model 

m mean 

ma mass averaged 

mer meridional 

min minimum 

n nozzle 

Pp plenum 

r radial direction 
ST stage 

t Ep 

u upstream of nozzle 
V vertical 

Z, axial direction 

9S tangential direction 
Superscriptss: 

a relative to rotor 
b exponent in deviation—angle relation 
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3. Design Specifications 
The fan is to deliver air at 8 = 1.3947 and 
R = 53.34 ft lb/lbm ” pes Design point conditions are as 
follows: 
a) Speed = 85,000 RPM 


b) Flow rate 83.5 lbm/min 


c) Driving power = 23.916 HP 
d) Conditions ahead of fan inlet-scroll=housing 
ares 
1) uniform flow 
2), ly ~waAnR 


3) P 13.656 psia 


Oo 


e) Required total pressure ratio between stator 
discharge and the entrance plane of the 
scroll-housing is P3/Py = 1.178 

The inlet-scroll—housing has a rectangular entrance 
(5.65 in. by 7.82 in.) and a flange 3.125 in. from the cen=- 
ter line of the fan. The distance of the center plane of 
the rotor from the back plate cannot exceed 2.1 in. because 
of the critical speed of the shaft. The hub diameter on 
the upstream side of the rotor cannot be less than 1.82 in. 
to accommodate the bearings. The discharge pipe of the 
fan has an ID of 3.75 inches. For manufacturing reasons, 
and to facilitate assembly, the outer diameters of the blad- 
ing should be constant. ‘The fan is to be single-stage and 


without inlet guide vanes. 
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4. Equipment and Instrumentation 

Air supply for the test facility is from a Carrier cen- 
trifugal compressor driven by a 300 HP electric motor. The 
compressor has a rated output of 4 lb/sec at a pressure 
ratio of 2. 

A general schematic of the test installation is shown 
in Figure 1. The major portion of the test equipment is 
shown in Figure 2. 

The duct from the plenum to the model is 4 ft. long and 
has 11.3 in. by 15.64 in. inside dimensions. It has two 
3 in. long pieces of honeycomb material to provide uniform 
flow at the entrance to the model. 

Flow rate through the model is controlled with the by- 
pass valves. The primary by<pass provides coarse adjustment 
of flow rate. The secondary by=pass, located in the test 
room, is used for fine flow rate adjustments. 

Flow rate is measured by a VDI standard flow nozzle, in 
the 4 in. inlet pipe, upstream of the plenum. The static 
pressure drop across the nozzle is measured on a 20 in. 
mercury micromanometer. Static pressure upstream of the 
nozzle is measured on a 92 in. water manometer. 

The plenum temperature and the temperature upstream of 
the nozzle are measured with iron—constantan thermocouples... 

Total pressure at the model entrance is measured by two 
Kiel probes mounted horizontally and vertically, and read 


on a 250 cm. water manometer. 
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Measurements in the model discharge annulus are made 
with a United Sensor three—dimensional five hole pressure 
probe. Pressure measurements are read On a 120 cm. water 
manometer. Yaw angle, YW , is read from an angle scale on 


the probe support. 
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5. Inlet Analysis 

Velocity distribution ahead of the rotor and pressure 
loss through the inlet must be known in order to start the 
fan design. This determination was made by testing a double- 
s@ale model of a proposed inlet design. The model casing is 
shown in Figures 3 and 4. The bellmouth for the model is 
shown in Figure 5. The belimouth hub is shown in Figures 6 
and 7. 

In order to generalize the test results, all flow rates 


were reduced to an equivalent flow rate 


“EE 


where p and T are the conditions at any point. 
Similarity Conditions 


The ability to apply data from the testing of one size 
machine to another is based on the concept of geometric and 
flow similarity. Geometric similarity was achieved by 
Changing all dimensions by the same ratio. Flow similarity 
is defined by the similarity parameters of Mach number and 
Reynolds number. The matching of these two parameters 
assures similar flow. 

Using the definitions of Mach number and the speed of 
sound, the Mach number, in terms of equivalent flow rate, 


becomes 


me ATS 


ES 


Noting that the area of the model is four times the area of 
the design and that % is the same, the equivalent flow 


rate of the model, for equal Mach numbers, is 
* * 
w (Model) = 4w (Design) 


The Reynolds number is defined by 


cae 
Nre = “ZT 


* 
and in terms of w it becomes 


i. i 


N,. = == 
Re ~ K VgrtT A 


Noting that the area of the model is four times the area of 
the design and that the characteristic length, 1, of the 
model is twice that of the design, the equivalent flow rate 


of the model, for equal Reynolds numbers, is 


Pp t 
w' (Model) = 2 fm —2 Mm w* (Design) 
Mp Pu tp 
Using the specification conditions for the design, and aver= 
age conditions for the model of 
Py = 14.7 psia 
550° R 


tm 
the equivalent flow rate of the model, for equal Reynolds 


numbers, becomes 
w" (Model) = 1.344 w" (Design) 


The equivalent flow rate for the design, calculated from 


the specifications, is 3.505 in.*. For Mach number simi- 


nbs 


larity this requires a flow rate of 6.8 lb/sec in the model. 
For Reynolds number similarity a model flow rate of 2.29 
lb/sec is required. 

Since the flow through the inlet is in the low subsonic 


range, the effects of Mach number are small and were neg= 
lected. The model testing was done at a flow rate which 
gave the Reynolds number similarity condition. 


Test Procedure 


- The flow rate was set and held constant by maintaining a 
constant static pressure upstream of the flow.measuring noz= 
Zle. Readings were taken at intervals during each traverse 


to determine the average flow rate and total pressure at the 


model entrance. 
The five hole pressure probe was traversed across the 
annulus height at the model discharge. Measurements were 


taken at each 10% point of the annulus height from the outer 
radius. Additional measurements were taken at:3, 6, 94, and 
97 percent of the annulus height to get an indication of -the 
boundary layer formation. Data was obtained to determine 


total and static pressure, and yaw and pitch angle. MTra- 


verses were taken at eight stations, spaced every 45 degrees 


around the annulus, as shown in Figure &. 
Data Reduction 


The reduced data for a typical model test is given in 
Appendix A. Due to the volume of data obtained only that 
for Modification 7 (Mod. 7) is given. 
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Equivalent flow rate at the nozzle was determined by the 
method presented by Kelly lila ° The relevant equation 
from {1 ] is 
wi = *€E ££ £ A VE AP. (5i1) 


n Py 


Equivalent flow rate at the model discharge is related 


to the equivalent flow rate at the nozzle by 


p Te 
* * Fy d 
WwW mm W orton \f nm be 


where pi, is the average static pressure at the discharge. 
Since the discharge velocities were small, the discharge 
static temperature was assumed equal to the total tempera- 
ture, which in turn was considered to be the measured ple= 
num temperature. 

The average static pressures at each of the eight sta- 
tions were determined from the pressure probe data. The 
data from the probe traverses was used in conjunction with 
the probe calibration chart (Fig. 9) to determine the 
following local flow conditions: the total and static 
pressures, yaw and pitch angles, and the local total= 


pressure—loss coefficient. The parameter used to enter 


Pa-Be 
Figure 9 1S) 6. - The further procedures are indicated 
a 1 . 2 
in Appendix A. 
The yaw angle W was measured directly from a vernier 
scale protractor mounted on the probe support. Positive yaw 


was defined as that producing clockwise whirl if looking 


lee 


upstream into the discharge annulus. The pitch angle 
was determined from Figure 9 and defined as positive for 
radial velocity component pointing outward from the dis- 
Charge centerline. 

In order to give an indication of loss characteristics 
for various points, stations, and model modifications, the 
local total-pressure=-loss coefficient was defined as 

7 al (5.3) 
where Py, is the average total pressure at the entrance 
to the model. 

In order to compare flow distributions at the various 
stations, and for different modifications, the dynamic 
pressure ratio a7a=-, was determined at each data point, 
where qa is defined as the theoretical dynamic pressure 
that should occur at the discharge annulus for uniform flow 


and a given discharge equivalent flow rate. That is, 


q=Pe-p (5.4) 
and 
ke 
ceed \ Bay (Sw 
dav ~ Ag 2 


Derivation of equation (5.5) was based on the following 


relations: 
continuity equation: w= 0... DAg Vay (5.6) 
equation of state: e == __fav__ (5.7) 
. av — ' 
g R eee 
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-.38) 


discharge equivalent k W i: \fR 
“ER ee US 


£Llow rate: Wa = 


definition of dynamic pressure: . Wavy 2 (5.9) 
(for incompressible flow) . av ° 


Substitution of V,,, from equation (5.6) into 


equation (5.9) gives 
eve ° 2 
_fav ; w 
Tav ~ “3° 012 g2 Age 20. 


av g2 A? 


or, from equation (5.7) 
°2 


Ke Rtsy 
Z 








g ss oe 
a 2g Ag Pav 


Substituting w from equation (5.8) gives 


= (wa") ———-— 
R 
Cay 2c Ag Pay 


. 2. 
=, eae 
day = A 


The dynamic pressure ratio then becomes 


therefore 


P oa 
= a Gey 
day (2. Pay 
= 2 


The local velocity parameter Viex/ Vay was considered use=- 





ful for noting the effectiveness of the various modifications 


made to the model in an attempt to attain axisymmetric: flow. 


zo 


For a local velocity V with components in the radial, 


axial, and tangential directions the axial velocity is 


V.,=Vcos © cos W . Therefore, for constant density 
Vax a 
2X \K.-- cos © cos YW (5.11) 
eae day 
Modifications 


Several modifications were made to the model in an 
attempt to attain the same flow conditions at each of the 
eight stations shown in Figure 8. Initial testing indicated 
that deficient flow existed at stations l, 2, and 8 Sta- 
tions 4, 5, and 6 had excessive flow, while stations 3 and 7 
were receiving the desired amount of flow. Yaw angles were 
excessive at stations 2, 3, 7, and 8 (about 115° 220m Pitch 
angles were highest at stations 4, 5, and 6 (about 89-10%). 

In an attempt to correct the flow pattern and achieve 
axisymmetric flow in the model discharge, seven internal 
modifications were made to the model. These modifications 
were tested individually except as noted below. The follow- 
ing is a description of each modification with a discussion 
of the reason for the change and a summary of the results 
obtained. For clarity, the pictures of each modification 
(except Mod. 7) are shown in the relevant figures without 


the hub and bellmouth assembled. 
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Mod. 1 -—- A flow splitter with a 2 in. radius of curva 
ture was placed at the top of the model casing (station 1) 
as shown in Figure 10. This change was made for two reasons. 
The original splitter was unavoidably damaged during the 
process of boring a probe hole for stations 1 and 5. It was 
also considered appropriate for the splitter to have some 
Curvature so as to provide a smoother change of flow direc-= 
tion at the top of the casing. The flow distribution 
(Fig. 11) remained unchanged and the local loss distribution 
(Fig. 12) was almost identical to the preliminary test re—- 
sults. For these reasons and due to the fact that this 
modification was used in all subsequent testing, the Mod. 1 
results were used as a basis for comparison with other 
modifications. 

Mod. 2 —- A flow splitter was placed at the bottom of 
the casing (station 5) as shown in Figure 13. It extended 
from the front wall of the casing to the beginning of the 
curvature of the bellmouth. This change was made to provide 
a smoother change of flow path and to direct more flow to 
the top of the casing. The only effect seen in this test 
was an increase in the pressure loss. The flow distribution 
was not altered. Therefore, these results are not shown. 

Mod. 3 = A full width flat guide vane was placed in the 
model entrance at the end of the flow splitter of Mod. 2. 
The guide vane alone is shown in Figure 14. It was set at 


zero angle of attack to the flow. The trailing edge has a 


circular cutout so that it wraps around a portion of the 
bellmouth. The maximum chord is 4.25 in. at the wall. 


PA 


This guide vane was designed and placed so as to block 
a portion of the flow going to stations 4, 5, and 6 and to 
force it to the upper stations 1, 2, and 8, where the flow 
was deficient. It was considered to be a gross over= 
correction as it was placed very close to the bellmouth 
lip and extended completely across the model entrance width. 
That is, an over=correction was made in hopes of later re=- 
ducing the correction to some optimum. 

The results for this modification showed that high pitch 
angles and high losses existed at station 5. The flow was 
deficient at this station near the outer radius but actual- 
ly increased in magnitude near the hub (inner radius). The 
flow conditions at station 1 were unchanged from those of 
Mod. 1. Therefore, these results are not shown. 

Mod. 4 — Two guide vanes were placed at the model en= 
trance as shown in Figure 15. The vanes are flat and have 
ai.5 in. chord. They were set at zero angle to the flow . 
and 2 in. from the plane of symmetry of the casing. This 
modification was designed to split the flow into three sec= 
tions. The middle section (between the vanes) would go to 
the lower stations 4, 5, and 6, and the two outer sections 
would go to the side and upper stations. The chord length 
was selected so as not to interfere with the flow close to’ 
the bellmouth. 

The supports for these vanes were designed such that each 
vane could be adjusted over a 2 in. slotted—track and set at 
any angle of attack. For this modification a neutral posi-= 


tion was used as a preliminary choice of vane positioning. 
24 


The results for this test showed no change in flow distribu- 
tion or losses from that of Mod. 1 and are therefore not 
presented. 

Mod. 5 ~ The guide vanes used in Mod. 4 were reset at a 
25 degree angle to the flow with the leading edge 0.6 in. 
from the plane of symmetry as shown in Figure 16. The re- 
positioning of the guide vanes in this manner was considered 
to be a large correction in a further attempt to reduce the 
flow going to the lower stations and to direct it to the 
upper stations. 

For the purpose of reducing experimental testing time 
such that other modifications could be tested, preliminary 
data for this modification was taken only for five points 
(10, 30, 50, 70, and 90 percent of the annulus height) at 
eaenm OL Stations 1, 3, 5, and sam if the results forecnr- 
data indicated significant improvement then a complete 
annulus traverse would have been carried out. 

Comparison to Mod. 1 can be made at the same percentage 
values of the annulus height. Results for Mod. 5 are shown 
in Figures 17 and 18 for axial velocity and loss distribu- 
tions, respectively. The flow at station 1 was increased 
slightly. The flow at station 5 was more uniform but of 
about the same amount as before. Changes at station 3 and 
7 were small. Figure 18 shows that the local losses in- 

. creased at station 5 and changed only slightly at stations l, 


3, and 7. 
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Mod. 6 = Two guide vanes were placed at the model en- 
trance as shown in Figure 19. The vanes are "S" shaped with 
a 5.35 in. chord. The leading edge was set at zero angle 
to the flow and was 0.75 in. from the plane of symmetry. 
From the leading edge they have a radius of curvature of 2 in. 
for 2.25 in., then a straight portion for 0.8 in., then an 
opposite radius of curvature of 3.5 in. for 2.55 in. to the 
trailing edge. The vanes are tangent to the bellmouth at 
the trailing edge. 

This modification was a gross correction to divide the 
flow into three sections and force sufficient flow to the 
upper stations. Although considered impractical this change 
was made and tested in the hope that an optimum and more 
practical configuration could be established. 

The Mod. 6 Velocity Distribution (Fig. 20) shows that in 
comparison with Mod. 1 the flow was increased at all sta- 
tions except 4 and 6 which were very deficient. Stations l, 
2, and 8 now had excessive flow, and the flow at station 5 
became even more excessive. Several stations had less uni- 
form axial velocity distributions than in Mod. 1. Stations 4 
and 6 experienced very high losses as shown in the Mod. 6 
Loss Distribution (Figt::21}4,%° 

Mod. 7 — A groove 1 in. deep and 2.28 in. long was cut 


in the bellmouth outer wall as shown in Figures 5 and 22. 


The purpose of this change was to provide more area for flow 


to the upper stations and to improve the bellmouth shape. 


The assembled configuration is shown in Figure 23. 
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Figures 24 and 25 show the axial velocity and loss dis- 
tributions, respectively. The results in general are very 
similar to those for Mod. 1. However, the station that 
Changed the most was station 1, and this was as an improved 
(increased) flow distribution. Stations 2 and 8 showed 
increased losses near the outer radius. Station 1 experi- 
enced somewhat lower and more uniform losses. The other 
stations indicated little or no change over Mod. l. The 
pitch and yaw angles were reduced slightly, but fairly sig- 
nificant yaw angles (15°) still existed at stations 2, 3, 7, 


and &. 


Predicted Conditions 


Ahead of Fan Rotor 


The final part of the inlet analysis was the prediction 
of the flow conditions ahead of the rotor for the fan design. 
This prediction was based on the non-dimensional velocity 
ratio and the total-pressure=-loss coefficient from the model 
tests, and on the fan=-inlet entrance conditions given in the 
specifications. 

Final Model Analysis. - The model test data indicated 
that the Mod. 1 and Mod. 7 configurations were similar. For 
better comparison the data for these two configurations was 


Circumferentially averaged for oS, VW ’ 4 , and Vee 
It was decided that Mod..7 gave slightly better flow distri- 


bution between the outer and inner walls, and that the 
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averaged yaw and pitch angles were somewhat smaller. The 
loss characteristics of the two configurations were nearly 
identical. Therefore, the Mod. 7 data was chosen for use 
in the prediction of rotor entrance conditions. 

The circumferentially averaged velocity distribution, as 
shown in Figure 26, was somewhat in error as the area under 
this curve was about three percent greater than unity 


(1.0316), although theoretically this area must equal unity; 


that is, 
4.0 
Vi... 
2% (z ) aio (5.12) 
O 
: av 
ag 
Vo 


Conditions that were considered likely causes of this dis- 
crepancy ure: decreased discharge area due to probe block= 
age, general instrumentation errors, improper estimation of 
the velocity profile near the walls, increased probe in- 
accuracies in the boundary layer, and the assumptions made 
in the calculations (such as tg = Tg , e = constant, etc.). 
TO correct this apparent error the Mod. 7 velocity dies 
tribution was normalized such that equation (5.12) was satis- 
fied. This was done by dividing each ordinate value of 
Vax/Vay by the area under the curve. Both distributions, 
Mod. 7 and corrected Mod. 7, are shown in Figure 26. 
The mass-averaged total-pressure-loss coefficient for the 


model was determined as follows 


4 7 Pr - Pma (5.13) 
MA Pma 7 Pma 


Ze 


where PL is the average total pressure at the model en- 
trance, and Py and Pha are the mass~-averaged total and 


static discharge pressures, respectively. tThus, 


sremcaliaiaie 
ma f aw 


which may be written as 


17 
eo! 
[ J es —= v., dr dg P 
us 
O 
Oo yg V 
av 
o Wk, 


For assumed constant Tre 


~T Ko 
[ J. r——= cdr d@pP 
ma 21T 
| i co ‘dr a@ 
9 


The static pressure p at discharge was mass-averaged 
in the same manner. The integration was carried out radi- 
ally and then circumferentially. Using equation (5.13) it 
was determined that 9). = 0.0484. 

Design Prediction. —- The conditions ahead of the rotor 
were determined by the method of Vavra [2 ] for an adia- 
batic expansion process with friction and varying cross- 
Section. 

The fan-inlet loss coefficient was assumed to be equal 
to the mass-averaged value for the model. 

The dimensionless flow function b was determined from 


the relation in [ 2 ] 


Zo 


(i eee . (5.14) 


where Ww, 1 and P, are given in the design specifica- 


0 
tions. The area A, ahead of the rotor is dependent upon 
the design choice since dimensional restrictions are given 
in the specifications. The double-sized model was based on 
the minimum hub distance and a casing diameter of 2.5 inches. 
Thus, the value of A, was prescribed. 

With ® and °? as known parameters, Table Cl of [2] 


was used to find the pressure ratio Po/P, . Thdés procedure 


was based on the simplifying assumptions that 


ye == 1.4 


and 
2b = 0.05 
The. following equations from [2] were then used to 


determine the velocity Vy : 


| Y=-1 
mike p . mM 
xX = = G 
O 
2 _ xlt% un X 
1l1- xX 


Xp 


Ss gi Rie 
| i 7a (8 ~1) nN oil 1/2 


mae [2- (2) ™] 


where Ny is defined as the polytropic exponent, Pp as the 


= 


polytropic loss coefficient, and X as a pressure relation 


used as a simplification of the equations. 
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The velocity V, represents the required average velo- 
city ahead of the rotor for a non-uniform velocity distri- 
bution, or the actual velocity for uniform flow that is 
Ora Oo le 


Since T) = ce 


for an adiabatic process, the static 
temperature was determined by 
2 


1 
1 l 
29 Ji Cp 


V 


The total pressure P, was determined by 


Egat ecw 
O 
Po 


I 
ale 


where 


ee oie 
2-(2 v1 


The design weight flow was used as a check on these pre- 


dicted conditions by 





ee A, V 
"D yay 1 ‘1 


The absolute Mach number is 
+5 
My) = Vi/(BSGR ty) 


Given below is a summary of the properties at the rotor 


Dole t s 
O 
P, = 13.534 psia T, = 714 R Vi = 627.6 ft/sec 
: O me 
py = 11.483 psia t, = 681.55" R ae 0.490 


a 
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Aerodynamic Fan Design 


The basic equations that are applicable to the design of 


an axial-flow compressor or fan are developed in [ 3]. 


These are formulated from the laws of conservation of matter, 


momentum, 


and energy, along with the thermodynamic equation 


of state of the working fluid. 


The simplifying assumptions that are made in [3] are 


as follows: 


are s 


ib) 


2) 


3) 


2) 


The general flow equations are applied only to 
compute flow distributions between blade--rows 
where blade forces are nonexistent. 

The flow is assumed to be steady and axially 
symmetric. 

Viscosity effects between blade rows are neg-= 
lected in the flow equations. However, these 
effects are considered by use of experimentally 
determined blade-element performance data. 
Empirical corrections are also made for wall 
boundary-layer effects on required flow area. 


Heat transfer is neglected. 


The two principal phases of the aerodynamic fan design 


1) 


2) 


Determination of blade=row velocity-diagrams 
for design point operation. 
Selection of blading to satisfy these velocity- 


diagrams. 
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Determination of Velocity-—Diagrams 
General. —- The determination of the various air veloci- 
ties and flow angles from hub to tip at the inlet and outlet 
of the rotor and stator was made based on the radial equili- 
brium equation [3] , which is 


2 
V dV QV 

Jgc RF = agt $5 + ve M+ 2+ ee L ; 

IP dr 3* Ox Sor iG ‘2 Se Vz Ox CS .) 


This equation results from the preceding simplifying assump- 
tions, and was the theoretical starting point for the fan 
design. 


It was assumed that the meridional velocity 


ae \ v,2 + vee is given by the axial velocity V, , 


and that YV-/dz = 0. Such a condition has been success- 
fully used in low=-aspect-ratio and lightly loaded blade=- 
row designs. Therefore, this assumption was considered 
appropriate for the fan design. 

It was decided that the design should include constant 
energy addition through the rotor blade-row. That is, the 
Change of total enthalpy through the blade-row is constant 
along a streamline. If the total enthalpy is assumed con- 
stant ahead of the rotor (Hy = constant), then a constant 


total enthalpy will exist behind the rotor. This may be 


(), « 


2 


expressed as 


The Euler turbine equation may be written as 
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a | [ (rv oe (rVe )4] (6.2) 
Since Vg = 0 ahead of the rotor, and Hyg - H, is constant, 
then (rVg ) 5 = constant. Thus, the above conditions re- 
guire that free-vortex flow exists behind the rotor blade- 
row. 
Therefore, neglecting radial velocity components, equa- 


tion (6.1) becomes 





2 
Ve Ve OVz 9S 
Ve “3s Pe Ve Jt te Ce 0 (6.3) 


The first two terms in this equation may be written as 


ve (ov 
-= (« — + ve ) 





or 
Vo ) (r Ve) 
"Yr ox 
Since rVo = constant, the sum of the first two terms in 


equation (6.3) equals zero. 


The radial equilibrium equation then reduces to 


Vv 
ve Se eR 0 (6.4) 


The simplest solution of equation (6.4) is the "simple- 
radial-equilibrium solution," based on the further assump- 
tion that the derivative of S with respect to r is zero. 
This condition has been widely used in compressor design, 
and was considered valid for the single-stage fan for slight 
change in annular area. This assumption is made only ata 


fixed value of @. 
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Therefore, equation (6.4) becomes 





with the solution 


dVz 
V = 0 
2G 
Vz = constant (6.5) 


Using this result as a basis, the velocity diagrams and 


flow properties were calculated for two assumed rotor inlet 


conditions; 


a) 


b) 


namely, 


Design A for the circumferentially-—averaged 
axial velocity for the most satisfactory flow 
determined from the inlet model tests. The 
masS—averaged total-pressure=-loss coefficient 
for the model was assumed. Total temperature 
and total pressure were assumed to be constant 
ahead of the rotor. Pitch and yaw components 
of velocity were assumed to be zero. 

Design B for uniform flow with axial velocity 
components only. The velocity and properties 
ahead of the rotor were based on the same loss 


coefficient as Design A. 


General calculation procedure. = The calculation method 


is one of iteration and compromise. Iteration is necessary 


in order that the continuity equation may be satisfied by 


proper annular dimensions. Compromise is necessary in order 


that blade-loading and Mach number limits are not exceeded. 


Experimental tests of axial-flow compressors show that a 


Mach number of 1.2 is the approximate limit for a transonic 
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eciagmemetnogd as USed in this paper. Further, the following 
diffusion—factor limits are appropriate for inlet stages: 
for the rotor tip, D less than 0.4; for the rotor hub, D less 
than 0.6; for the stator, D less than 0.6 [ 3 ] - Although 
these criteria are somewhat vague, they remain the best 
indicators of design limitations that are currently avail- 
ebilc 

The detailed development of the calculation method is 
presented in [3] . ‘The iterative procedure is started by 
selecting the radial stations according to percentages of the 
passage depth. The computation scheme and results for velo- 
city calculations are shown it Tables I and II, for Design A 
and Design B, respectively. Only the final results and the 
final calculation lineups of the iteration procedure are 
presented. 

The velocity leaving the stator was assumed to be entire- 
ly axial, thus providing maximum flow rate for a given stage 


loading. 


Design A. = Table I shows the steps taken to determine 
each parameter. Figure 27 shows a general stage velocity- 
diagram which is indicative of the fan design. 

From the results of the inlet-model analysis, the assumed 


inlet conditions to the rotor were as follows: 





T, = 714 R 

Pi = 13.534 psia 

V 

ie as shown in Figure 26 
av 
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iat 627.6 ft/sec 


sd 
II 


lt 1.74 in. (for an assumed tip clearance of 0.010 
: in.) 

Tih: Oo gen 

The conditions at station 1 are listed as steps 1 through 
8 in Table I. The iteration was started by assuming that 
the hub radius was constant through the rotor (step 9). Most 
of the steps are clearly indicated under the procedure col- 
umn in Table I, in which Figure 27 was used. However, sev-= 
eral steps had to be related to important design assumptions 
and will be discussed in detail. 

Yo 2 (step 11) was determined by use of the Euler rela- 


tion, equation (6.2), and the definition of DP, the driving 


power. That is, 


DP = w(H5 - H)) 
Wey 
or DP SS eV ea) 
Jg rea 


Since w = TTN/30, the tangential velocity behind the 
rotor may be found as a function of fro. 

It has been pointed out in [3] that large axial velocity 
reductions across any blade-element are not conducive to high 
efficiency. In addition, the available experimental aaiee. for 
transonic rotors shows that the range of Vz,2 to Ve1 ratios 
across the rotor are as follows: for Ehe tip, 0.85 to 1.05; 
for the mean radius, 0.90 to 1.15; for the hub, 0-95 to Wsav. 
Therefore, in step 14 the axial velocity through the rotor 
at the mean radius was assumed to be constant. This assump- 
tion causes the axial velocity changes near the tip and hub 


regions to be typical of the experimental data. 


er 


Steps 24 through 26 involved a compromise in which sever- 
al important design conditions were considered. These were 
as follows: the critical blade=-chord Reynolds number, fric- 
tional losses, shock losses, blade-loading, manufacturing 
restrictions, and the practicality of rotor diameter measure=- 
Te Waters 

The critical blade-chord Reynolds number may be defined 
as the limiting value of Reynolds number below which a large 
rise in loss occurs. Above this value the losses vary only 
Slightly or not at all. The available experimental data are 
generally for flow conditions at or aboveithis critical 
Reynolds number. Therefore, rotor blade design that is 
based largely on test data should also be for similar flow 
conditions. Based on the relative inlet velocity to a 
blade-element, Vavra [4] indicates that this limiting 
Reynolds number is about 2 x 10°. 

Too high a solidity is associated with increased fric- 
tional losses. Too low a solidity results in higher blade 
loading, and increased shock losses near the tip. Recom- 
mended rotor tip solidity for transonic rotors in regard ame 
shock losses is 1.0 [3 ]. 

For manufacturing reasons the minimum blade spacing at 
the hub is about 0.25 inches. 

An even number of rotor blades allows for easier measure- 


ment of the tip diameter. 


Therefore, in consideration of the above facts, the mini- 


5 


mum blade-chord Reynolds number was set at 2 x 10° for the 


a6 


rotor. The corresponding minimum blade-chord was deter- 


mined, for the known conditions ahead of the rotor, by 


VV, *€ 
Nee = aw &e (6.6) 
TA 


where 4, was found in [ 5 ) for known values of t, - 
several combinations of chord distributions and blade 
numbers were investigated in order to determine a satisfac-— 


tory solidity range. Pertinent equations used were 


n|a 


It was considered important to have a tip solidity fairly 
close to 1.0 as recommended in [ 3 ] , but also essential to 
establish a solidity range similar to that for the experi- 
mental data for transonic rotors. 

As a result, 12 rotor blades with uniformly varying 
chord were selected such that the Reynolds number was at or 
above 2 x nile A summary of this analysis, which was based 


On no area change through the rotor, is as follows: 


‘ a ’ 


¢ 


Radial position Vy Cada os S or 
TIP 1450.1 OmooS OST 0 OvoET 0.823 
..LO% 139525 O'. 557 0.760 0.868 0.876 
50% 169s 0.665 0.800 0.694 ass 
90% 960.4 0.809 0.840 Oma20 1 eS 
HUB 914.0 0.850 0.850 0.476 1.786 


From this summary the linear chord distribution was ex- 


pressed as a function of radius for those iterations to fol- 


low. Thus, 
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It should be noted that the above chord-solidity summary 
refers only to the first iteration. With reduced blade 
height at the rotor discharge the values of Yr, and Yo5 
will be slightly different for a given percentage of passage 
depth. A streamline will follow a curved path through the 
rotor, depending une annulus constriction. Therefore, the 
solidity must be averaged for the values at rotor entrance 
and exit. This was done in step 27, where a, is referred 
to the chord and spacing at station 1 for a given percentage 
passage depth from the casing, and T> to the chord and 
spacing at station 2 for the same percentage passage depth. 
The total-pressure-loss coefficient wo was determined 


from the experimental data for a geometrically similar 


double=circular-arc transonic rotor fe] - This rotor was 
é 


tested under similar flow conditions. The values of W 
were for minimum=-loss incidence angle. This method of 
determining the loss coefficient was considered appropriate 
for the simplified theory used in this paper. 

The relations used in steps 30 to 32 are developed in 
Sl: 

It was considered probable that the boundary layer in the 
fan rotor would have significant blockage effects, especially 
due to the small blade heights. Therefore, in step 34 a 
blockage factor of 0.96 was assumed. The continuity equation 


was checked by graphical integration. 
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The iteration was repeated for a new hub radius at sta# 
tion 2, reestarting at step 9 until continuity was satisfied. 
The final rotor outlet conditions were used as stator 

inlet conditions. A similar process was followed for the 
stator design. However, the stator blade solidity was deter- 
mined in a somewhat different manner than for the rotor, as 
excessive solidities resulted for the criterion that minimum 
Reynolds number be equal to 2 x 10>. 

A solidity distribution was sought which would be simi- 
lar to the available experimental data for the NACA~-65 (A, 9) - 
series airfoil. (The specific reasons for choosing this 
airfoil shape will be given later). The extensive test data 
for this airfoil includes solidities up to about 1.6, but 
extrapolated data for blade solidities up to 2.0 are given 
met ol . 

Several combinations of blade-chord and number of blades 
were investigated. As a result, the solidity at the 90% 
position (near the hub) was set at 1.55, and the solidity 
at the 10% position (near the tip) was set at 1.0. Blade 
chord was chosen as constant at 0.805 inches, and the number 
of blades was set at 13. 

Step 42 indicates the resulting Reynolds numbers. These. 
values were deemed acceptable for two reasons [3] . First, 
the relative difference between Reynolds numbers of 1.5 x 10° 


ana 2.0 x 10° 


is small for the NACA-65 (Aq) - series airfoil 
under moderate blade loadings. Secondly, lightly loaded 


blades (such as the fan design of this paper) have only a 
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gradual increase in loss below the critical Reynolds number, 
rather than the more typical steep rise in loss for normal 

blade loadings. The diffusion factor values listed in step 
43 indicate that the fan stator is rather lightly loaded as 
compared to the suggested design limit for stators of D less 


than 0.6 . 


The total-pressure-loss parameter (step 44) was deter- 
mined from [3] for known diffusion factors. 

The continuity equation was checked using an assumed 
hilockage factor of 0.94. Iteration was found to be unneces-= 
sary. 

The stage total pressure ratios and stage adiabatic effi- 
ciencies for the three radial positions were determined. 
These values were mass averaged as shown in Table I. 

Finally, the total pressure ratio over the inlet-fan 
combination was determined for comparison with the design 


specification. 


Design B. = Table II shows the procedural steps taken 
to determine each parameter. 
From the results of the inlet analysis for uniform flow, 


the conditions ahead of the rotor were taken as: 
Vi = Vl of DESIGN A = 627.6 £t/sec (uniform) 


t, = 681.59 R 
Total properties and rotor radii the same as Design A. 
The detailed steps and design assumptions were similar to 
Design A. Iteration and mass-averaging were likewise similar- 


ly done. 
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selection of Blading 


General. —- The blade-sections at each of the three 
radial positions were selected in order to satisfy the de- 
sign velocity-diagrams. 

Double=-circular-arc blading was chosen for the rotor be- 
cause this relatively simple airfoil shape has been used 
effectively in transonic compressors. 

The NACA=65 (A, ))-series airfoil was considered to be 
appropriate for the stator blading. The reasons for this 
choice were the magnitudes of the design inlet Mach numbers 
for the stator and the availability of extensive experiment- 
ali-Gata for thats airfoil. 

The parameters that describe blade-element flow are 
total-pressure loss, incidence angle, and deviation angle. 
The total-—pressure-loss coefficient for minimum loss inci- 
dence has already been established in the determination of 
the velocity~diagrams. Therefore, the design incidence 
angle is for minimum loss. 

Camber angle, incidence angle, and deviation angle are 
used to define the blade camber, air approach, and air 
leaving directions, respectively. These angles are shown in 


Fig. 28 for a typical blade-element. 


Rotor. —- A rotor=-blade selection procedure is presented 
i) [3] for the double-circular-arc airfoil. This method 
is based on a correlation of experimental test results with 


both annular and rectilinear cascades. Empirical correction 
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factors are used to modify the basic two-dimensional-cascade 
predictions for design blade shape. These factors are based 
largely on the testing of single-stage transonic compressors. 
Therefore, this method was considered appropriate for the 
fan design. The following equations from [3] were used 


to determine the blade=shape parameters: 


AB’ =P +O-S- (6.8) 
a Ki(<,),. + ng (Cees) (6.9) 
f= K5(So). + p +(e Co te $) +/( + (5 “S.-o) (6.10) 
D = nB'+ (SeSi2p)-(ée&0)]I -($¢ $5) l-k; (ishotKs(Se),6 a 
I-™m™ +7 
—t 


The additional unknown parameters on the right side of 
these equations were determined from the velocity-diagram 
properties and an assumed blade thickness distribution as 
follows: for the rotor tip, t/c equal to 0.050; for the ro- 
tor mean radius, t/c equal to 0.075; for the rotor hub, t/c 
equal to 0.100. The thinning of the blade toward the tipe 
was considered appropriate for the supersonic velocities pre- 
sent there, and the associated shock losses. In addition, 
this particular distribution matches that for the transonic 
rotor which was used for the determination of the radial 


loss distribution [6]. 
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The procedural steps for the determination of camber, 
incidence, and deviation angles are shown in Tables III and 
IV for Design A and Design B, respectively. 

For both designs the calculated camber angle for the 10% 
radial position was found to be negative. This occurred due 
to the very low turning angle required. For this case, a 
zero-camber blade section was chosen at this location. The 
corresponding values of incidence angle and deviation angle 
were determined by equations (6.9) and (6.8), respectively. 
Only the corrected values are shown in Tables III and IV; 


namely, those for a Zero camber angle. 


Stator. -—- The stator blade-sections were determined 
from the two-dimensional-cascade performance data for the 
NACA-65 (A,,)-series airfoil as presented in [3]. a 
blade selection method based on a correlation between ob- 
served compressor performance and cascade predictions was 
not followed, as available stator data were considered to 
be too limited to establish reliable correlations. 

The following equations from [3] were used to deter- 


mine camber, incidence, and deviation angles: 


AB - [ (Co), 2 (Sey) 


—— (6.12) 
[1 -me+n] - - 


i a Gp ~ Cr + nYPp (6.13) 
° 2 = mM 
fo = So = (SK). + 7% P (6-44 
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The additional unknown parameters on the right side of 
these equations were determined from the velocity-diagram 
properties and an assumed thickness distribution of t/c 
equal to 0.10 for the entire blade height. The choice of 
10% maximum thickness was considered appropriate for the 
stator velocity-diagram conditions. This choice also simpli- 
fied the use of the cascade performance curves in [3] . 

The procedural steps for the determination of camber, 
incidence, and deviation angles are shown in Tables III and 


IV for Design A and Design B, respectively. 
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7. Results and Discussion 


Inlet Model 


The primary result of the inlet model testing was the 
fact that axisymmetric flow could not be achieved for modifi- 
cations that would be practical and have reasonably¥ low 
losses. 

The various modifications have been discussed in sec- 
tion 5. The axial velocity distributions for Mod. 1, Mod. 5, 
Mod. 6, and Mod. 7 are shown in Figures 11, 17, 20, and 24, 
respectively. The loss distributions for Mod. 1, Mod. 5, 
Mod. 6, and Mod. 7 are shown in Figures 12, 18, 21, and 25, 
respectively. Results for the other modifications were con- 
sidered insignificant and are not shown. 

Mod. 7 was selected for the fan design because it gave 
the best flow distribution of the configurations tested. 
However, the improvement over the original inlet-model was 
Slight. The Mod. 7 results for axial velocity distribution 
were circumferentially averaged as shown in Figure 26, and 
corrected for system errors such as probe blockage and in- 
strumentation inaccuracies. This averaging process was 
necessary so that available fan design procedures could be 
used. However, Figure 26 could present a misleading impres- 
sion of the inlet-model test results and should be appraised 
in conjunction with Figure 24 which shows the Mod. 7 distri- 


bution by stations. 
Most of the modifications produced only small changes of 


the flow pattern. Only the large "S" vanes of Mod. 6 


- 
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produced large variations from the original flow distribu- 
tion. This modification caused excessive flow to exist at 
the previously deficient stations 1, 2, and 8. The losses 
for these stations were also lowest for Mod. 6 as compared to 
all other configurations. It seems possible that a modifica- 
tion similar to Mod. 6 could produce axisymmetric flow, but 
only at the expense of high losses. Such a configuration 
could evolve from a trial and error process of reducing the 
surface area, curvature, leading-edge spacing between vanes, 
and/or chord of the "S" vanes. 

The reason for the lack of axial symmetry for more prac- 
tical modifications such as Mod. 7 is largely due to the 
area constriction» forced on the flow as it approaches Mime 
top of the model casing (near station 1). The slopes of the 
Casing (13 degrees) as shown in Figure 4 are excessive. 
Another reason for the lack of axial symmetry might be that 
the top of the model was not extended far enough from the 
annulus centerline. This distance was not specified but 
would necessarily be a design consideration for an aircraft 
system. 

The curvatures of the model bellmouth contours are ex- 
cessive. This was indicated by the thicker boundary-layers 
and higher velocities near the outer radius for most of the 
eight stations. The larger velocities form an adverse gra- 
dient ahead of the fan; that is, the relative inlet Mach 


number at the rotor tip is increased. This condition is es- 
pecially critical as this Mach number is an important design 
Limi te 
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The Mod. 7 pitch and yaw angles are shown with the data 
reduction procedures in Appendix A. Pitch angles were great- 
est at stations 4, 5, and 6 where they averaged about ~9° 
(radially inward). Yaw angles were highest at stations 2, 3, 
7, and 8 where they averaged about 17 degrees. 

It should be noted that the algebraic signs of the pitch 
angles as shown in Appendix A are taken with respect to the 
pressure probe. Positive pitch angle is for flow that is 
radially outward. Thus, the signs for © for stations l, 
2, 7, and 8 are correct as indicated; and the signs for 
stations 3, 4, 5, and 6 must be reversed to comply with the 
established sign convention. Yaw angles are shown with 
proper signs. That is, positive yaw angles occur for tan- 
gential velocity components producing clockwise flow if look- 


ing upstream at the discharge. 
Fan Design 


The results for Design A and Design B are very similar 
because the averaged velocity distribution for Design A is 
nearly uniform (Fig. 26). The performance figures and 
blade-section values listed below are for both designs. 

The velocity-diagram determinations for Design A and 
Design B are shown in Tables I and II, respectively. The 
selection of blade-sections are shown in Tables III and IV 
for the designs in the same order. 

The adiabatic stage efficiency (based on total conditions) 


was 90.5 percent, and the total-pressure ratio was 1.243. 


49 


The latter value exceeds the design specification require~ 
Teme Onmwieivs Dy >.> percent. 

Twelve double-circular-arc blades were selected for the 
rotor. The blade chord varied uniformly from 0.85 in. to 
0.75 in. from hub to tip, respectively. The maximum thick- 
ness distribution varied uniformly from 0.10 to 0.05 from 
hub to tip, respectively. The blade height at rotor en- 
trance was 0.83 inches. The blade height at rotor discharge 
was 0.74 inches. 

Thirteen NACA-65 (A, 9) -series blades were selected for 
the stator. A constant blade-chord of 0.805 in. and a con- 
stant maximum thickness ratio of 0.10 were used. The height 
of the stator blades was constant, equal to 0.74 inches. 

The rotor blade camber angle was zero degrees near the 
tip (10% location). This occurred due to the low air turn- 
ing angle in that region (A sa 

It was necessary to circumferentially average the axial 
velocity distribution for Design A as the design method is 
based on the assumption of axial symmetry. The results of 
this averaging gave nearly uniform flow as shown in Figure 
26. However, Figure 24 shows that the actual flow distribu- 
tions are lacking in axial symmetry. Therefore, the Design A 
results may be misleading. In addition, the large changes 
Of flow angle for Mod. 7 would cause large variations from 
the design incidence angle for the rotor blades. This con- 


dition is especially undesirable in transonic rotors as low 
losses occur only over a very narrow range of inlet flow 


angles at high Mach numbers as pointed out in [3] : 
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However, Robbins and Glaser [7] have shown that cir- 
cumferential inlet flow distortion of velocity and flow angle 
may not have as significant an effect upon compressor per- 
formance as might be expected. In [7] the results of 
considerable circumferential distortion showed that the adi- 
abatic efficiency and pressure ratio were unchanged, and 
that the equivalent weight flow was slightly reduced. These 
results pertain to tests on a five-stage transonic compres- 
sor. After the first stage most of the flow angle distortion 
became equalized and the velocity distortion was reduced by 
about 30 percent. 

Therefore, the lack of axial symmetry in the Design A 
inlet may not be critical. However, the compact size of the 
fan and the fact that it is only single stage makes com- 
parison with [ 7 less significant. 

Thus, it seems apparent that actual testing would be 
necessary to determine the effect of the inlet-distortion on 
the predicted fan performance of Design A. The Design B 
blading was designed for uniform flow and could be used for 
a comparison if tested with the Design B assumed rotor en- 
trance conditions (or ‘similar’ flow conditions). 

The désign results show that the curvature changes through 
the rotor blading were small. That is, the three streamlines 
that were considered (10%, 50%, and 90% locations) experi- 
enced relatively small radius changes. Therefore, the design 
assumption of zero entropy gradient behind the rotor was con- 


sidered valid. 
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The simplifying assumptions that V, equals Vmer and 
dv,/) z = 9 were also considered appropriate because the 
design blading had low aspect ratio (about 1.0) and was 
lightly loaded, having small diffusion factors. 

A critical area in the determination of the velocity dia- 
grams was the chord-solidity selection, as the flow through 
the blade=-rows was close to the critical blade-chord- 
Reynolds number. This condition occurred due to the com- 
pactness of the fan. 

Operation at higher blade-chord-—-Reynolds number could be 
attained by re-design. Since the diffusion factors are 
relatively low the number of rotor and stator blades could 
be reduced by one. For the same solidity level the chord 
could be increased in proportion to the blade spacing in- 
crease. Thus the blade-chord-Reynolds number would be higher 
and the flow would be less likely to experience losses 
associated with low Reynolds number. 

This paper presents a simplified approach for the aero- 
dynamic fan design which is commonly used for the design of 
axial-flow compressors. The complicated real flow effects 
such as secondary flow, tip clearance losses, and unsteady 
flow were ignored, as they were considered to be beyond the 
scope of this study. Furthermore, it is realized that for 
an actual design it might be necessary to repeat the design 
calculations to assure a satisfactory se aise that is sound 


both from an aerodynamic and a mechanical point of view. 
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Table III. Blade-Sections for Design A 


(a.) Rotor 
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Parameter Procedure 
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Table III. (Concluded) 
(b.) Stator 


Radial Position 
Parameter Procedure 
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Table IV. Blade-Sections for Design B 


(a.) Rotor 
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Parameter Procedure 
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